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Résumé 


La géostatistique est utilisée pr étudier les variations spatiales de propriétés fonctionnelles du transfert de l'eau, en 
relation avec l'organisation des sols. Une approche expérimentale est menée sur un bassin versant de 10 ha de Bretagne 
dont les sols sont développés dans des limons sur schiste briovérien. 


Dans un premier temps, la teneur en eau pondérale, dont la mesure est facile et précise, est la variable utilisée pour 
aborder des aspects méthologiques : influence de l'échelle de l'échantillonnage et du contexte pédologique sur la forme du 
variogramme. 


Dans un second temps, des cartographies par krigeage de la teneur en eau pondérale sont analysées en relation avec 
la carte des sols, présentant notamment les variations spatiales de l'hydromorphie des sols. Cette démarche analyse la 
signification fonctionnelle des limites pédologiques et la variabilité spatiale de la propriété étudiée à l'échelle des unités 
pédologiques. Elle concourre a l'étude hydrologique du bassin versant. 


Enfin le transfert de l'eau in situ est décrit en différents sites d'une parcelle, à l'aide du traçage isotopique. La 
géostatistique permet, dans cette troisième application de quantifier des paramètres du transfert de l'eau, d'en apprécier la 
variabilité spatiale et de mettre en évidence des hétérogénéités spatiales que la carte des sols n'avait pas laissé prévoir. 


Abstract 


Geostatistics has been used to study the spatial variations of the functional properties of water transfer in relation to 
soil characteristics. An experimental approach has been applied to a 10-hectare drainage basin in Brittany where the soils 
are derived from loëss lying on Brioverian schist. 


The weight percent water content, is easily and accurately measured and is the variable used in a first step to deal 
with the methodological aspects, i.e. the influence of the sampling scale and soil characteristics on the variogram. 


Secondly, maps of the gravimetric water content established by kriging are analyzed in comparison With the soil map, 
showing in particular the spatial variations of soil hydromorphy. In this step the functional significance of the soil 
boundaries and the spatial variability of the property studied, at soil unit level are analyzed. This assists in the hydrologic 


study ofthe drainage basin. 


Lastly, in situ water transfer is described at different places on a plot of land, using isotope tracing. Geostatistics here 
enables the quantification of water transfer parameters, an assessment of their spatial variability and the detection of 


spatial heterogeneities not suggested by the soil map. 


Introduction 


L'étude du fonctionnement hydrique des sols in 
situ est classiquement abordée par deux types 
d'approches : dans la première, une approche 
pédologique, on se sert des caractéristiques 
morphologiques du sol pour déduire, de façon 
qualitative, le fonctionnement hydrodynamique du 
sol dans l'espace ; dans la seconde, une approche 

hysique, on étudie quantitativement ce 
onctionnement sur quelques stations à l'aide de la 
mesure de différents paramètres et de modèles 
physiques. L'utilisation de la géostatistique permet 
de concilier et de valoriser mutuellement ces 2 
approches, prenant en compte les aspects spatiaux et 
quantitatifs du fonctionnement hydrodynamique. 
Son application, développée dans différents domaines 


des sciences de la terre [J.P. Delhomme, 1976], est 
récente en science du sol [T.M. Burgess et R. Webster, 
1980 ; E. Byers et D.B. Stephens, 1983 ; Y.M. Gajem 
etal., 1981 ; A.B. McBratney et R. Webster, 1983 ; M. 
Vauclin et al., 1983] ; elle est adaptée aux propriétés 
hydriques du sol qui sont des variables régionalisées, 
c'est-à-dire présentant un caractère à la fois aléatoire 
et structuré dans l'espace. La géostatistique permet, 
d'une part, d'analyser la variabilité spatiale d'une 
propriété donnée, en particulier à l'aide du 
variogramme, fonction rendant compte de la 
covariance entre les points de mesure. Elle permet, 
d'autre part, de reconstituer au mieux la réalité, à 
partir d'un échantillonnage et d'apprécier la fiabilité 
des estimations faites par krigeage. 
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Fig. 1. - Variogrammes des données de teneur en eau mesurée sur un 
transect à 20 (0), 40 (A), 60 (-) et 80 cm t°), niveaux correspondant 
respectivement à des horizons pédologiques différents 


L'approche choisie ici comporte deux étapes : la 
première porte sur les modalités de description de la 
variabilité spatiale des propriétés hydriques du sol 
dont l'importance a été soulignée par D.R. Nielsen 
[1973] ; la seconde analyse les variations spatiales de 
quelques propriétés hydrodynamiques, en liaison 
avec la pédologie, rendant compte de l'organisation 
des sols dans l'espace, l'hydrologie et le fonc- 
tionnement hydrodynamique d'un terrain 
expérimental. Différents aspects de l'utilisation de la 
géostatistique en physique du sol sont ainsi illustrés 
ici. 


Terrain expérimental 


Le cadre expérimental choisi sur un bassin 
versant élémentaire de 10 ha, situé à 30 km au sud- 
est de Rennes -à l'ouest de la France-, le bassin 
versant de Nouvoitou. Il s'agit d'une zone bocagère, 
peu pentue ; les sols sont développés dans des limons 
sur schiste briovérien. Des études préalables menées 
sur ce site, ont permis une connaissance de la 
distribution spatiale des sols par une cartographie à 
1/1 000 et du fonctionnement hydrique des sols sur 
quelques stations ayant fait l'objet d'un suivi tensio- 
neutronique. L'étude a porté sur la période hivernale 
durant laquelle le fonctionnement hydrique des sols 
est varié dans l'espace et peut être mis en relation 
avec l'organisation de sols différenciés latéralement 
et de façon privilégiée par des critères d'hydromor- 
phie. Ainsi sur les sols sains de versant, l'infiltration 
de l'eau se fait verticalement ; vers l'aval, où se 
développe une hydromorphie de surface, il existe une 
composante latérale de sub-surface. 


Résultats 
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Fig. 2. - Variogrammes des données de teneur en eau mesurée à 20 


cm de profondeur, à différentes échelles 


1 - Analyse de la structure d'une variable 
donnée 


La structure d'une variable donnée est étudiée à 
l'aide du variogramme, fonction rendant compte de 
comment se détériore l'information acquise en un 
point, lorsque l'on s'en éloigne. L'interprétation est 
généralement faite pour la variable étudiée ; 
cependant, on a cherché à préciser si l'échelle et le 
contexte pédologique influençaient les résultats. 


La variable étudiée est la teneur en eau 
pondérable, variable dont la mesure est simple et 
précise. Ce paramètre, relativement stable dans le 
temps sur le versant, peut être assimilé à la capacité 
au champ ; vers l'aval, on passe de façon différenciée 
suivant les niveaux à une quasi-saturation. 
Différents plans d'échantillonnage ont été envisagés : 
(1) un transect de 3 m situé à l'aval du bassin versant ; 
(2) un transect de 30 m situé sur le versant ; (3) une 
parcelle de 2,5 ha dont les points de l'unité de versant 
ont été uniquement retenus ; (4) une parcelle de 2,5 
ha dont tous les points de mesures ont été retenus ; (5) 
le bassin versant entier. La teneur en eau a été 
mesurée à 20, 40, 60 et 80 cm. Les variogrammes ont 
été construits pour chacun de ces 5 plans 
d'échantillonnages et de ces 4 niveaux [C. Gascuel- 
Odoux, 1984]; seuls certains sont ici présentés. 


L'influence du contexte pédologique est illustré à 
la figure 1, où sont comparés les variogrammes 
obtenus à 20, 40, 60 et 80 cm à partir de l'échan- 
tillonnage du transect de versant situé sur un sol 
brun faiblement lessivé : ces 4 niveaux correspondent 
à des horizons de natures pédologiques différentes. 
Ainsi le variogramme est pépitique à 40 cm, comporte 
une portée de l'ordre de 4-5 m, à 20 et 60 cm, et 
présente une allure linéaire à 80 cm. Des résultats 
analogues sont obtenus à partir des autres échan- 
tillonnages. Au-delà même de la variable étudiée, il 
est donc nécessaire de prendre en.compte de façon 
étroite le contexte pédologique de la mesure : unité 
pédologique, horizon. 


L'influence de l'échelle est illustrée à la figure 2 
où sont comparés les variogrammes obtenus à partir 
de la mesure de la teneur en eau à 20 cm sur les 
différents échantillonnages. Ceci confirme l'existence 


mêtres 
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Fig. 3. - Cartographie par krigeage de la teneur en eau, mesurée à 
40 cm de profondeur ; comparaison à la position de la limite A/A2, 
t—) 


d'effets d'échelle, soulignée par ailleurs par Y.M. 
Gajem et al. [1981] et P.A. Burrough [1983]. 
Cependant, il semble que la variabilité spatiale soit 
identique avec des structures légèrement différentes 
si on se limite à l'unité pédologique. Il existe donc des 
effets d'échelle allant à l'encontre d'une géné- 
ralisation des résultats obtenus à une échelle donnée. 


2 - Analyse de caractéristiques de la 
couverture pédologique 


Dans l'étude du fonctionnement hydrique d'un 
terrain expérimental, la géostatistique permet de 
décrire les variations spatiales d'une propriété 
donnée et de les interpréter en regard du tracé des 
limites pédologiques. Cette approche précise donc le 
caractère fonctionnel des limites pédologiques et 
valorise une étude ponctuelle. Sur le bassin versant 
étudié, on a cherché à relier les variations de la 
teneur en eau pondérale, mesurée sur les différents 
échantillonnages décrits ci-dessus, aux limites pédo- 
logiques d'hydromorphie, limites dépendant de 
l'importance et de la durée de l'engorgement des sols. 


Ainsi la cartographie de la teneur en eau mesurée 
à 40 cm de profondeur à l'échelle de la parcelle ou du 
bassin versant, permet de caractériser la limite 
A2/A9g, représentant l'apparition de l'hydromorphie 
au sein de l'horizon lessivé, Av (fig. 3) : à l'amont de 
cette limite, la teneur en eau apparaît constante 
malgré une situation géomorphologique de versant ; 
la limite coïncide avec une augmentation rapide et 
importante de la teneur en eau ; en aval de cette 
limite, la teneur en eau croît vers le centre du bassin 
versant. D'autres analyses faites à différentes dates 
confirment ces résultats, montrant une liaison étroite 
et stable au cours de la saison hivernale entre les 
caractéristiques d'hydromorphie et la teneur en eau. 


Par contre, au sein des horizons de surface et des 
horizons de profondeur, les limites pédologiques sem- 
blent nettement moins corrélées aux variations de la 
teneur en eau. Ainsi, à 20 cm de profondeur, la teneur 
en eau présente, au cours de la saison hivernale, des 
variations temporelles qui expliquent les difficultés à 
établir des relations entre la morphologie des 
horizons et leur teneur en eau. Comparant le com- 
portement des niveaux 20 et 40 cm, on peut donc dire, 


qu'il existe, au sein du niveau de surface, un écou- 
lement latéral de sub-surface. 


Ainsi des propriétés fonctionnelles du transfert de 
l'eau et des caractéristiques morphologiques du sol 
repérées lors de la cartographie sont analysées 
conjointement, non pas à l'échelle de quelques 
stations, mais à celle de l'unité ou de la limite 
pédologique. Cette approche permet, d'une part, 
d'analyser la signification fonctionnelle de traits 
pédologiques, en terme de transfert ; d'autre part, d'en 
tirer les conséquences au niveau de l'hydrologie du 
bassin versant étudié. Une des limites de cette 
approche est de ne prendre en compte qu'une variable 
donnée. 


3 - Analyse du transfert de l'eau à 
l'échelle d'une parcelle 


Le transfert de l'eau dans le sol est abordé, in situ, 
à l'aide du traçage isotopique. L'oxygène 18 est le 
traceur utilisé, traceur météorique intrinsèque, c'est- 
à-dire provenant des précipitations et marquant la 
molécule d'eau elle-même. En effet, les variations de 
la concentration en oxygène 18 des précipitations 
dépendent des épisodes de condensation et d'éva- 
poration que l'eau subit : ainsi les précipitations 
hivernales, de septembre à mars, constituent un 
signal isotopique classiquement utilisé pour suivre la 
propagation de l'eau dans le sol [B. Blavoux, 1978 : 
R.O. Munnich, 1983]. 


Cette technique est particulièrement adaptée à 
l'étude des variations spatiales du transfert : le signal 
d'entrée est identique à l'échelle de la parcelle et 
facilement reconstitué à partir de prélèvements 
réguliers des précipitations (fig. 4a) ; les réponses du 
sol à ce signal sont reconstituées deux mois et demi 
après le pic négatif de janvier, et en différents sites 
sur lesquels des carottages permettent de recons- 
tituer les profils isotopiques (fig. 4b). 


L'analyse géostatistique a porté sur deux 
paramètres du traçage, rendant compte des modalités 
du transfert de l'eau : la profondeur et l'amplitude du 
pic des profils isotopiques qui, comparés au signal 
d'entrée, rendent respectivement compte de la vitesse 
effective du transfert de l'eau et de la dispersivité du 
milieu. Sur l'unité de sol brun faiblement lessivé 
étudiée, on met en évidence une vitesse moyenne 
effective de 8,7.10-8 m.s-1 ; le coefficient de variation 
est de 30 % et le variogramme est pépitique : les 
variations spatiales sont donc moyennes, mais 
rapides dans l'espace. Par contre, l'analyse des 
données d'amplitude du signal, ayant trait à la 
dispersivité du sol, montre une organisation spatiale 
de la parcelle : il existe un net gradient que la carte 
des sols n'avait pas laissé prévoir (fig. 5). 


Cette application illustre l'intérêt de la 
géostatistique pour analyser la variabilité spatiale 
d'une fonction de transfert donnée. On peut ainsi 
obtenir la valeur moyenne sur la parcelle, le 
coefficient de variation et mettre en évidence des 
hétérogénéités spatiales. Le traçage isotopique, 
permettant de décrire le transfert effectif de l'eau 
dans le sol, est un outil analytique avantageux 
lorsqu'il s'agit d'étudier le transfert de l'eau en de 
nombreux sites. 
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Fig. 4. - Signal d'entrée du traçage (A) : concentration en oxygène 18 et hauteur d'eau des précipitations hivernales. Réponses du sol (B) : 


profils isotopiques, situés dans le plan d'échantillonnage, au 12 mars 


AVAL 


Fig. 5. - Cartographie par krigeage de l'amplitude du signal isotopi- 
que, exprimé en 0%ol80 et rendant compte de la dispersivité du milieu, 
sur une parcelle de sol cartographiécomme homogène 


Conclusion 


Les propriétés hydriques du sol présentent dans 
l'espace un aspect chaotique tel qu'il est difficile 
d'analyser l'information acquise en de nombreux 
points. La géostatistique, utilisée ici dans différents 
objectifs, s'avère un outil intéressant pour décrire la 
variabilité spatiale d'une propriété donnée, pour 
intégrer les résultats à une échelle donnée et pour 
donner une représentation cartographique des 
hétérogénéités spatiales. En science du sol, cette 
ps P est en particulier valorisée en regard de la 
carte pédologique, donnant une description spatiale 
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de l'organisation des sols et servant de point d'appui à 
l'étude du fonctionnement hydrique. Les résultats 
obtenus sur un terrain expérimental, concernant la 
description de la variabilité spatiale, semblent 
dépendre à la fois de l'échelle de l'échantillonnage et 
du contexte pédologique. Une généralisation de cette 
approche du fonctionnement hydrique passe par une 
prise en compte des effets d'échelle, une meilleure 
connaissance de l'organisation des sols et un élargis- 
sement à l'étude conjointe de‘plusieurs variables. 
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